Transport in a One-Dimensional Superfluid: Quantum Nucleation of Phase Slips 
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We present an analytical derivation for the quantum decay rate of the superflow through a weak 
link in a one-dimensional (ID) Bose-Einstein-condensate. The effective action for the phase dif- 
ference across the link reduces to that of a massive particle with damping subject to a periodic 
potential. We find an algebraic flow-pressure relation, characteristic for quantum nucleation of 
phase slips in the link and show how short-wave length fluctuations renormalizing the interaction 
between Bosons remove the quantum phase transition expected in this class of systems. 
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Quantum fluids exhibit fascinating phenomena as they 
push quantum mechanical effects to the macroscopic 
level. Classic examples are the electronic condensate 
in superconductors and the uncharged bosonic ''He and 
fermionic '^He superfluids Furthermore, the new 

laser- and evaporative cooling techniques have opened up 
new opportunities in the fabrication and manipulation of 
Bose-Einstein condensates (BEC) in various alkali-metal 
vapors . A spectacular phenomenon is the friction free 
transport, often studied in terms of the Josephson effect 
across a constriction or orifice separating two reservoirs. 
Here, we concentrate on a specific geometry, a narrow 
channel, and study the decay of the superflow due to 
quantum generation of phase slips. Such decay processes 
have been studied in superconducting wires, both exper- 
imentally 1^ and theoretically In this letter, we in- 
vestigate the quantum decay of driven uncharged con- 
densates — the Galilean invariance and the Schrodinger 
dynamics in these systems pose new challenges in the the- 
oretical description and produce drastic changes in the 
physical results as compared to the charged situation. 

In a (thick) one-dimensional superconducting wire, 
topological fluctuations of the phase (phase slips) are 
weak, resulting in a true superconducting response ^] 
with an algebraic current-voltage (I-V) characteristic 
V oc J^ff~^ and a dimensionless coupling constant g w 
SOrg/Ai (rg denotes the radius of the wire and Xl the 
London penetration length). Reducing g, e.g., via de- 
creasing the diameter of the wire, the system undergoes 
a (T = 0) quantum phase transition at 5 = 2 to a metal- 
lic phase with an ohmic characteristic. In the metallic 
wire the dynamics of the phase derives from integration 
over fermionic modes and is determined through the cou- 
pling to the electromagnetic field. Here, we study an un- 
charged BEC which is accurately described by the Gross- 
Pitaevskii Lagrangian for the condensate wave function 
ip{x,t). The Galilean invariance of the homogeneous sys- 
tem prohibits the nucleation of phase slips and the homo- 
geneous system is always in the superfluid state. Includ- 



ing a weak link breaks Galilean invariance and phase slips 
can nucleate at the defect ||^, however, we find that the 
link always remains in the superfluid state and the entire 
system exhibits no quantum phase transition whatsoever. 
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FIG. 1. The ID channel is connected to two superfluid 
reservoirs with a weak link at a:: = 0. The superflow / in 
the channel is driven by the phase difference A<1? = — "l>_t 
between the reservoirs. 

We start with the Gross-Pitaevskii Lagrangian giving 
the dynamics of the condensate wave function 1). With 
m the mass of the bosons, po the condensate density, and 
U the strength of the repulsion, the Lagrangian reads 



u 



Cgp = ifi tpdrip - i^(dxi'){,dxil}) - ^ - Po] ■ (1) 

The system admits states ijj — ^/pj ex-p{imvx/'h) car- 
rying a flow / = piv with the renormalized density 
Pi = Po ~ Tnv'^/2U and the velocity v, driven by the 
phase difference A$ = mvL/h across the reservoirs. 
The suppression of the modulus produces the critical 
flow Ic = {2 / 3)'^^'^ pq y^U Po / m. A crucial feature of the 
Lagrangian (^ is its Galilean invariance: with ipi^^ t) 
a state, further solutions are generated by the boost 
exp(iTO [—vx + /Tijiplx — vt, t), with v the relative 

velocity of the two coordinate systems. A flow carrying 
state is mapped via such a boost to a stable state with 
7 = 0, hence homogeneously flowing states are stable. 
An inhomogeneity in the tube breaks the symmetry and 
quantum nucleation of phase slips bound to the perturba- 
tion appear. Here, we introduce a weak link (see Fig. |^) 
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modeled by a variation in the stiffness, 



(2) 



where c{x) vanishes for |a;| > a/2 (a denotes the width 
of the weak hnk). The action of the channel then reads 



S^Sr + Sl 



dt 



dx [Cgp + Cw] 



(3) 



with S'fl(L) the contribution of the right (left) lead. Note, 
that the above perturbation is time independent with the 
weak link at rest with respect to the reservoirs. This 
contrasts with Ref. where the flow around a mov- 
ing object has been considered. Applying the mapping 
ip{x,t) e'}cp{—iinvx /?i)ip{x — vt,t) and pi allows 

for a comparison of our system with the setup in 0| . 

In order to study the quantum behavior of the super- 
fluid we go over to the Euclidean action H 



Se = 



in/3 



dr i dx\fiil)drip + 'H\}l','il^'\ 



(4) 



with T the imaginary time, f3 = 1/T the inverse temper- 
ature, and Ti. -0] the Hamiltonian density, 



A saddle point solution extremizing the above action in- 
volves ^ and Tp as independent variables; this freedom 
allows for an analytic continuation of the integration con- 
tour in the path integral. Performing functional deriva- 
tives produces the equations 

-hdrip^ -{n^/2m)dltp + [/ [ v^v - po] V' , 

hdr^=~{h^/2m)dl'^+U[^^p-po]^. (5) 

Note that with ^P{t), '0(r) satisfying Eq. (||) the config- 
uration tp *(— t), iP*{—t) is also a solution, with ip* the 
complex conjugate of tp. Furthermore, the periodicity of 
"0 and ^{x,t -\- hj3) = r), and the boundary 

condition ?l/3/2) = ^p{x, fif3/2) imply that 



(6) 



for the saddle point solution. 

The ID quantum phase slips, saddle point solutions 
describing the quantum tunneling of the ^ field, take the 
topological form of vortex-antivortex pairs in the (x, t)- 
plane [Q. Within the vortex core region the modulus of 
the wave function drops to zero, while in the region far 
from the vortex center its variation is small and the ac- 
tion is determined by the behavior of the phase. In the 
uncharged superfluid, Galilean invariance prohibits for- 
mation of a vortex-antivortex pair in the homogeneous 
region and thus the vortices are bound to the weak link. 



In the following we derive an effective action for the phase 
difference across the link. 

We first integrate out the weak link. For a small width 
a < ^ the dominant contribution in the action is the ki- 
netic term and the equation for the wave function ip 
inside the weak link simplifies to 



a,{[i-c(x)]a,V^} = o 



(7) 



The boundary conditions are 'p{a/2) ~ ')jjj^(a/2) and 
ip{—a/2,T) — ipL{~a/2) with ipn [ipL] the wave func- 
tion in the right (left) lead. Inserting the solution of (0) 
into the action we obtain the contribution 
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~h 



dT-^\rpL'pL+1pR'4>R-TpR'>pL-^L'>pR\ ■ (8) 

^Po 



Note, that in the above expression tppi and -pL are taken 
aX X = a/2 and x — —a/2. The critical flow of the weak 
link I^^ = {m/fipo) /°^^2 '^■^ ~ c(a;)] ^ is assumed to 
be small compared to the critical flow Ic of the homoge- 
neous system I^/Ic ^ 1- 

Next, we focus on the contribution of the leads. The 
flow in the wire is limited by the small critical current Iw 
of the link, justifying the ansatz 



(9) 



with h, dxp, and small, of order Iw/Ic- For the fol- 
lowing calculations, it is convenient to define the length 
^ = %/ y/rnpoU , the sound velocity v1 = poU/m, and 
the time tq — £,/vs- Inserting the ansatz (||) into (||) , the 
modulus and phase decouple to lowest order in Iw/Ic, and 
using (^ the equation for the phase (p takes the form 

- To drP){x,T) + {e/2) dlP{x,T)^sin ,t)\ , (10) 

with P~{x,t) = [P{x,t) — (j){x,—T)]/2. The left side de- 
fines a diffusion equation, while the right side describes 
a nonlinear and nonlocal coupling of the wave function. 
The boundary conditions demand that the asymmetry 
vanishes in the limit a; — > oo and r ±?i/3/2, see 
Eq. (||). The solution of the linearized system then is a 
promising ansatz for describing the solution far from the 
vortex core region. Below we find that the asymmetry is 
of order Iw/Ic outside the core region and the relevant 
behavior of the phase slip is determined by the region 
where the linearization of Eq. (|lo| ) is justified. Apply- 
ing a Fourier transformation leads to a local and linear 
system of equations for p{uj) — u{lu) + iv^Lu) 
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(11) 



defining modes with a dispersion relation —uP'/v^ = 
A:^ (l + $^fc2/4). Next, we solve (0) for the phase in 
the right lead x > a/2 with the boundary conditions 
0fl(a^ cOi r) = and (f>R{a/2, uj) = pioj) + iqitLi), where 
V'fl(x,T) = exp{ipR{x,T)), 
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- [A+p + 2cc;rog]e-^-(2---)/2?| , 



v(x,uj) = — ; — < e 

^ ' ' A_ - A+ \ 2wToA 

_X+p + 2ujTo q 

2cjToA-|_ 



(13) 



-s_(2a;-a)/25 



with A± = 2(1 ± Y^l — w^Tq ) and s± = ^/\±. Inserting 
the solution back into the effective action Se, Eq- @), 
we concentrate on the first term (the second term in- 
volves the Hamiltonian H = J dxTL which is conserved 
under time evolution and therefore does not contribute). 
Linearizing in the small quantities h, dx(t>^ and 9r'/>, the 
variation h of the modulus drops out and the action is 
determined by the phase field alone, 

/°° dio 
— / dx Lo u{XtUu) v{XtUu) . (14) 
-oo Ja/2 

A straightforward calculation leads to the effective action 
for the phase <j)j^ {a/2, r) on the right side of the weak link 



Se 



hK 



duj 



TT 7o 47r (1 + WTo)i/2 



1 + UJTo 



2q^ 



with K = TTpaS^ the number of particles per coherence 
length. An expansion in ujtq is equivalent to an expan- 
sion in Iw/Ic and we obtain to lowest order 



Se — / ctw [ IwIp — 2ujpq\ 

^ J — oo 



(15) 



The first term is a Caldeira-Leggett type damping pX| ] 
for the symmetric part p of the phase (j)ji{a/2,T) and 
describes the fiow of energy from the weak link to the 
reservoirs via sound wave excitations. Such a term is ex- 
pected for a dynamics admitting sound waves and is also 
present in the charged superconductor. The second term 
is a coupling between symmetric and asymmetric parts 
and is only present for Schrodinger dynamics. 

An analogous calculation with 0l(— a/2, t) adds a sim- 
ilar effective action for the phase on the left side of the 
link. Adding up terms from the leads and the link we 
obtain the total effective action for the phase difference 
(/j(t) = <j)ii{a/2,T) — 0i(— a/2,T) across the weak link 
(minimization with respect to the 'center of mass' vari- 
able z9(t) = (j)B{a/2,T) + (j>L{-a/2,T) gives ■& = 0). In- 
troducing the notation = [(/3(t) ± (/?(— r)] /2 for the 
symmetric and asymmetric parts we obtain 



^ = ^ldr 

h 2-K 



dr 



,b+(r)-^+(r')]^ 



-if (t)9^(^+(t) 



An (r - t') 

+ /dT{/„ [l-^-(r)V2-cos</7+(T)] -I^+{t)}. (16) 



Integrating out (p~ in Scs produces the mass term 
h J dr mw{drf^Y /2 with m^, — K"^ /AI^-k^ via the re- 
lation if^ ~ —{K/27tI^) Collecting terms, the 
effective action for the symmetric part v?"*" of the phase 
difference across the link is equal to the action for a par- 
ticle with mass rriuj and damping 77 — K/2tt in a periodic 
potential and driven by the force /. The general form of a 
phase slip is ip~^{T) = (t) + g~ (t) + a.TCsin{I / 1^,) where 
the shift described by last term is due to the driving force 
/, while {g~) denotes a kink (anti-kink) associated 
with the presence of a vortex (anti-vortex) in the link. 
Applying the trial function = ±2 arctan [(2t ± r)/ Tc] 
with the kink width Tc — K/nlw [0, the action for the 
phase slip takes the form 



5cff /h = 2K In (t/tc) - 2ttIt - In ?/ , 



(17) 



where Iny is found by collecting terms independent of 
T. The kinks interact logarithmically with strength 2K, 
while the driving force / produces a linear repulsion 
within a pair. Using the above trial function g*, we 
can justify the linearization in (p^: the region of large 
asymmetry cj)~(x^T) is limited by ^ along x jl^ and by 
the width of the kink along the r direction. The con- 
tribution from the core then is independent of r and only 
enters through the fugacity y. 

Thermodynamics [I — Q): In the following we consider 
a gas of quantum phase slips. The associated action de- 
scribes charged particles with logarithmic interaction in 
ID and a chemical potential Iny. The one-dimensionality 
of the gas follows from the confinement of the vortices to 
the weak link. This contrasts with the superconducting 
wire where phase slips appear throughout the homoge- 
neous channel and the action of the quantum phase slip 
gas is mapped to that of 2D charged particles with loga- 
rithmic interaction of strength 2g Q . 

In the latter system, a quantum phase transition (at 
g = 2) has been found 0| as quantum fiuctuations drive 
the system ohmic at small values of g < 2 (including a 
weak link and ignoring phase slips in the wire shifts this 
transition to g = 1 [T^ ). We then may expect a simi- 
lar transition to occur in the superfiuid a,t K = 1. In- 
deed, a careful analysis shows that a massive particle with 
damping and subject to a periodic potential undergoes a 
quantum phase transition for a critical damping K = 1 
]l^ . However, in the above derivation of the effective ac- 
tion ( |l6| ) we have neglected high frequency fiuctuations 
of the order parameter which affect the long wave length 
behavior of the homogeneous superfiuid via a renormal- 
ization of the sound velocity Vs = \/ f {'j)pQh/m, where 
7 — mU/h^po- The universal function 7(7) is obtained 
from the exact (Bethc Ansatz) solution of interacting ID 
Bosons [p^. 
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As the original relation Trpoh/m = Kvg between the 
sound velocity and the dimensionless parameter K is not 
renormalized, a consequence of the Galilean invariance 
, the renormalization of K can be found from that of 
the sound velocity. For 7^1, the result K = TrpoC ^ 1 
is recovered, while in the limit 7 — > cxd we obtain hard- 
core bosons with K = I For intermediate values of 
7, K varies smoothly between these two limits. There- 
fore, no quantum phase transition is present in the un- 
charged superfluid. 

Superflow (/ 7^ 0): The shape of the phase slip is deter- 
mined by the driving current /: minimizing the (T — 0) 
action ( |l7| ) with respect to t we find the kink-anti-kink 



separation r = K/nl (see Ref. 11 1 for a detailed dis- 
cussion addressing the motion of a particle with damp- 
ing in a periodic potential). The kink width Tc deter- 
mines the asymmetry of the phase slip via the relation 
ip^ oc —drf'^, see Fig. |[ This form of the phase slip 
agrees well with numerical simulations 




FIG. 2. Phase slip solution: the dotted line is the symmet- 
ric part tp^ , while the dashed line represents the asymmetric 
part . The width of the asymmetry is determined by Tc 
and the distance between the kink - anti-kink pair is r. 

The quantum decay rate of the superflow is determined 
by r = B exp{~Sps/h) Q with Sps the action of the 
phase slip and B arising from fluctuations around the 
saddle point. At T = 0, the action is dominated by the 
logarithmic interaction between the kink and the anti- 
kink leading to a factor I^^ . Gaussian fluctuations in the 
distance between the kink and anti-kink with frequency 
oc / lead to a preexponential factor B (x l/I and there- 
fore F oc (in the superconducting wire the relative 
positions of the kink - anti-kink pair fluctuate both along 
X and T and therefore B cx 1//^, leading to F oc 
Note, that the prefactor produces a sizeable correction of 
the flow rate as compared to the saddle point result [|j . 
At finite temperature we have to determine the transition 
rates F* both to lower and higher energy states. These 
transitions then produce a time dependent phase differ- 
ence between the reservoirs which in turn determines the 
drop in the chemical potential across the channel, 

2K-1 



= 27r F+ -F" 



/ 

T 

hi,„ 



2K-2 ^ (19) 



With the pressure p oc A/Lt we arrive at an algebraic flow- 
pressure characteristic for T = 0, while at T ^ the 
channel always exhibits a linear behavior for a small su- 
perflow /. 

In the end, we have found a remarkable similarity in 
the response of charged and uncharged condensates: the 
quantum decay of the superflow due to phase slips leads 
to algebraic current-voltage and flow-pressure relations, 
respectively. On the other hand, we flnd distinct differ- 
ences as well: first, the quantum nucleation of phase slips 
in the present uncharged situation is bound to a defect, a 
consequence of Galilean invariance, while the same phe- 
nomenon in the charged case can take place throughout 
the wire. Second, the presence of high energy modes 
renormalizing the sound velocity removes the quantum 
phase transition for the uncharged superfluid, preserving 
the superflow. 

We wish to thank K. Le Hur and R. Heeb for valuable 
discussions. 
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